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ABSTRACT 

Our velocity relative to the rest frame of the cosmic microwave background (CMB) generates a dipole temperature anisotropy on the 
sky which has been well measured for more than 30 years, and has an accepted amplitude of v/c = 1.23 x 10"^, or v = 369 kms"'. 
In addition to this signal generated by Doppler boosting of the CMB monopole, our motion also modulates and aberrates the CMB 
temperature fluctuations (as well as every other source of radiation at cosmological distances). This is an order 10"^^ effect applied 
to fluctuations which are already one part in roughly 10"^, so it is quite small. Nevertheless, it becomes detectable with the all-sky 
coverage, high angular resolution, and low noise levels of the Planck satellite. Here we report a first measurement of this velocity 
signature using the aberration and modulation effects on the CMB temperature anisotropics, finding a component in the known dipole 
direction, {l,b) = (264°, 48°), of 384km s"' ± 78kms-' (stat.) ± llSkms"' (syst.). This is a significant confirmation of the expected 
velocity. 

Key words. Cosmology: observations - cosmic background radiation - Reference systems - Relativistic processes 



1. Introduction 

This paper, one of a set associated with the 2013 release of 
data from the Planck'' mission (Planck Collaboration I 2013), 
presents a study of Doppler boosting effects using the small scale 
temperature fluctuations of the Planck CMB maps. 



* 'And yet it moves," the phrase popularly attributed to Galileo 
Galilei after being forced to recant his view that the Earth goes around 
the Sun. 

** Corresponding author: D. Scott <dscott@phas . ubc . ca> 

t Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



Observations of the large cosmic microwave background 
(CMB) temperature dipole are usually taken to indicate that our 
Solar System barycentre is in motion with respect to the CMB 
frame (defined precisely below). Assuming that the observed 
temperature dipole is due entirely to Doppler boosting of the 
CMB monopole, one infers a velocity v - (369 ± 0.9) km s ' in 
the direction (l,b) = (263.99° ± 0.14°, 48.26° ± 0.03°), on the 
boundary of the constellations of Crater and Leo (Kogut et al. 
1993; Fixsen et al. 1996; Hinshaw et al. 2009). 

In addition to Doppler boosting of the CMB monopole, ve- 
locity effects also boost the order 10"^ primordial temperature 
fluctuations. There are two observable effects here, both at a 
level of y6 s v/c = 1.23 x 10"^. First, there is a Doppler "modu- 
lation" effect, which amplifies the apparent temperature fluctua- 
tions in the velocity direction. This is the same effect which con- 
verts a portion of the CMB monopole into the observed dipole. 
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The effect on the CMB fluctuations is to increase the amplitude 
of the power spectrum by approximately 0.25% in the velocity 
direction, and decrease it correspondingly in the anti-direction. 
Second, there is an "aberration" effect, in which the apparent ar- 
rival direction of CMB photons is pushed toward the velocity 
direction. This eff'ect is small, but non-negligible. The expected 
velocity induces a peak deflection of p - 4.2' and a root-mean- 
squared (rms) deflection over the sky of 3', comparable to the 
effects of gravitational lensing by large-scale structure, which 
are discussed in Planck Collaboration XVII (2013). The aber- 
ration effect squashes the anisotropy pattern on one side of the 
sky and stretches it on the other, effectively changing the angu- 
lar scale. Close to the velocity direction we expect that the power 
spectrum of the temperature anisotropies, C(, will be shifted so 
that, e.g., ^ = 1000 ^ 1001, while ^ = 1000 ^ ^ = 999 in the 
anti-direction. In Fig. 1 we plot an exaggerated illustration of the 
aberration and modulation effects. For completeness we should 
point out that there is a third effect, a quadiTipole of amplitude 
0^ induced by the dipole (see Kamionkowski & Knox 2003). 
However, extracting this signal would require extreme levels of 
precision for foreground modelling at the quadrupole scale, and 
we do not discuss it further. 

In this paper, we will present a measurement of /3, exploiting 
the distinctive statistical signatures of the abeiTation and mod- 
ulation effects on the high-^ CMB temperature anisotropies. In 
addition to our interest in making an independent measurement 
of the velocity signature, the effects which velocity generates on 
the CMB fluctuations provide a source of potential bias or con- 
fusion for several aspects of the Planck data. In particular, ve- 
locity effects couple to measurements of: primordial "tnl "-type 
non-Gaussianity, as discussed in Planck Collaboration XXIV 
(2013); statistical anisotropy of the primordial CMB fluctua- 
tions, as discussed in Planck Collaboration XXIII (2013); and 
gravitational lensing, as discussed in Planck Collaboration XVII 
(2013). There are also aspects of the Planck analysis for which 
velocity effects are believed to be negligible, but only if they are 
present at the expected level. One example is measurement of 
/NL-type non-Gaussianity, as discussed in Catena et al. (2013). 
Another example is power spectrum estimation; as discussed 
above, velocity effects change the angular scale of the acous- 
tic peaks in the CMB power spectrum. Averaged over the full 
sky this effect is strongly suppressed, as the expansion and con- 
traction of scales on opposing hemispheres cancel out. However 
the application of a sky mask breaks this cancellation to some 
extent, and can potentially be important for parameter estima- 
tion (Pereira et al. 2010; Catena & Notari 2012). For the 143 
and 217 GHz analysis mask used in the fiducial Planck CMB 
likelihood (Planck Collaboration XV 2013), the average lensing 
convergence field associated with the aberration effect (on the 
portion of the sky which is unmasked) has a value which is 13% 
of its peak value, corresponding to an expected average lensing 
convergence of /? x 0.13 = 1.5 x lO"'*. This will shift the angular 
scale of the acoustic peaks by the same fraction, which is de- 
generate with a change in the angular size of the sound horizon 
at last scattering, 0, (Buries & Rappaport 2006). A 1.5 x 10""* 
shift in 0» is just under 25% of the Planck uncertainty on this pa- 
rameter, as reported in Planck Collaboration XVI (2013) — small 
enough to be neglected, though not dramatically so. Therefore it 
does motivate us to test that the observed aberration signal is not 
significantly larger than expected. With such a confirmation in 
hand, a logical next step is to correct for these effects by a pro- 
cess of de-boosting the observed temperature Notari & Quartin 
(2012); Yoho et al. (2012). 




Fig. 1. Exaggerated iflustration of the Doppler aberration and 
modulation effects, in orthographic projection, for a velocity 
V - 260 000 km s"' = 0.85c (approximately 700 times larger 
than the expected magnitude) toward the northern pole (indi- 
cated by meridians in the upper half of each image on the left). 
The aberration component of the effect shifts the apparent posi- 
tion of fluctuations toward the velocity direction, while the mod- 
ulation component enhances the fluctuations in the velocity di- 
rection and suppresses them in the anti-velocity direction. 



Before proceeding to discuss the aberration and modulation 
effects in more detail, we note that in addition to the overall pe- 
culiar velocity of our Solar System with respect to the CMB, 
there is an additional time-dependent velocity effect from the or- 
bit of Planck (at L2, along with the Earth) about the Sun. This 
velocity has an average amplitude of approximately 30kms"', 
less than one-tenth the size of the primary velocity effect. The 
abeiTation component of the orbital velocity (more commonly 
referred to in astronomy as "stellar abeiTation") has a maximum 
amplitude of 21" and is corrected for in the satellite pointing. 
The modulation effect for the orbital velocity switches signs be- 
tween each 6-month survey, and so is suppressed when using 
multiple surveys to make maps (as we do here, with the nominal 
Planck maps, based on a little more than two surveys), and so 
we will not consider it further here. 

2. Aberration and modulation 

Here we will present a more quantitative description of the aber- 
ration and modulation effects described above. To begin, note 
that, by construction, the peculiar velocity, j8, measures the ve- 
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locity of our Solar System barycentre relative to a frame, called 
the CMB frame, in which the temperature dipole, a\„,, vanishes. 
However, in completely subtracting the dipole, this frame would 
not coincide with a suitably-defined average CMB frame, in 
which an observer would expect to see a dipole Ci ~ lO"'", 
given by the Sachs-Wolfe and integrated Sachs-Wolfe effects. 
The velocity difference between these two frames is, however, 
small, at the level of 1 % of our observed v. 

If T' and h ' are the CMB temperature and direction as 
viewed in the CMB frame, then the temperature in the observed 
frame is given by the Lorentz transformation (see, e.g., Challinor 
& van Leeuwen 2002; Sollom 2010), 



r(n) 



r(n') 



(1) 



(2) 



r(i-H p)' 

where the observed direction n is given by 

M ' -H [(y - l)n ' ■ V + yySJv 
r(l+n'-j8) ' 

and y = (\ - ff-y'^^^. Expanding to linear order in (3 gives 

T'in ') = T'(h - V(n ■ P)) = To + 6T'{n - V(m ■ fi)), (3) 

so that we can write the observed temperature fluctuations as 

6T{n ) = Ton ■ )8 + 6T'(n - V(/i ■ )8))(1 + n ■ /3). (4) 

Here To = (2.7255 + 0.0006) K is the CMB mean temperature 
(Fixsen 2009). The first term on the right-hand side of Eq. (4) is 
the temperature dipole. The remaining term represents the fluc- 
tuations, abeiTated by deflection V(n ■ fi) and modulated by the 
factor (1 + « ■ )8). 

The Planck detectors can be modelled as measuring differ- 
ential changes in the CMB intensity at frequency v given by 



Iv(v, n ) 



2hv' 



1 



exp[hv/kBT(n)] - I 



(5) 



We can expand the measured intensity diff'erence according to 



6Iv(v, n ) 



dly 
df 



6T{n) 



1 (f% 

2 dT^ 



6T^(n) + . 



(6) 



To 



Substituting Eq. (4) and dropping terms of order (ff) and (6T^), 
we find 



6Iy(v, h) - 



'df 



{Ton ■p + 5T'{h'){\+byn -P)), (7) 



where the frequency dependent boost factor by is given by 



by - — coth 



(2V0) ^' 



(8) 



with V() = ksTo/h - 57 GHz. Integrated over the Planck band- 
passes for the 143 and 217 GHz channels, these effective boost 
factors are given by = 1.961 + 0.015, and b2n = 3.071 ± 
0.018, where the uncertainties give the scatter between the indi- 
vidual detector bandpasses at each frequency. We will approxi- 
mate these boost factors simply as 143 = 2 and b2n - 3, which 
is sufficiently accurate for the precision of our measurement. 
The inferred temperature fluctuations will then be 



6Iy(v, n ) 
dIy/dT\ 



Ton ■p + 5T\n - V(n ■ )8))(1 + byh ■ yS). (9) 




Fig. 2. Decomposition of the dipole vector p in Galactic co- 
ordinates. The CMB dipole direction {l,b) = (263.99°, 48.26°) 
is given as jSy, while two directions orthogonal to it (and each 
other) are denoted as p^ and p^. The vector j8x lies within the 
Galactic plane. 



Notice that, compared with the actual fluctuations in Eq. (4), the 
modulation term in Eq. (9) has taken on a peculiar frequency de- 
pendence, represented by by. This has arisen due to the coupling 
between the fluctuations and the dipole. Ton ■ y8, which leads 
to a second-order term in the expansion of Eq. (6). Intuitively, 
the CMB temperature varies from one side of the sky to the 
other at the 3 mK level. Therefore so does the calibration factor 
dly/dT, as represented by the second derivative d^Iy/dT'. We 
note that such a frequency dependent modulation is not uniquely 
a velocity effect, but would have arisen in the presence of any 
sufficiently large temperature fluctuation. Of course if we mea- 
sured T(n ) directly (for example by measuring Iy{v, n ) at a large 
number of frequencies), we would measure the true fluctuations, 
Eq. (4) (i.e., we would have a boost factor of by = 1). 

3. Methodology 

The statistical properties of the aberration-induced stretching 
and compression of the CMB anisotropics are manifest in "sta- 
tistically anisotropic" contributions to the covariance matrix of 
the CMB, which we can use to reconstruct the velocity vec- 
tor (Buries & Rappaport 2006; Kosowsky & Kahniashvili 201 1; 
Amendola et al. 2011). To discuss these it will be useful to in- 
troduce the harmonic transform of the pecuHar velocity vector, 
given by 



dnYl^{n)p- 



(10) 



Here (3im is only non-zero for dipole modes (with L - 1). 
Although most of our equations will be written in terms of Plm, 
for simplicity of interpretation we will present results in a spe- 
cific choice of basis for the three dipole modes of orthonormal 
unit vectors, labelled jSy (along the expected velocity direction), 
)8x (perpendicular to jSy and parallel to the Galactic plane, near 
its centre), and the remaining perpendicular mode p^. The direc- 
tions associated with these modes are plotted in Fig. 2. 

In the statistics of the CMB fluctuations, peculiar velocity 
effects manifest themselves as a set of off-diagonal contributions 
to the CMB covariance matrix: 



{TfimiT (21112) CMB - y^/~l)^ ( 
LM ^ 



(x h L 
nil m2 M 



To 



(2A + 1)(2^2 + 1)(2L+1) „ 

X ^/ — Wf,t2L/^LM, (11) 
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where the weight function W^'' is composed of two parts, related 
to the aberration and modulation effects, respectively: 



(12) 



The aberration term is identical to that found when considering 
gravitational lensing of the CMB by large scale structure (Lewis 
& Challinor 2006; Planck Collaboration XVII 2013), 



(1(21 



I + ^_lYt+t2+L 



j ^JL(L+lyl(^l + l) 



xC 



,TT / <^ 1 



£2 L 
1 -1 



+ (A^^2), (13) 



while the modulation term is identical to that produced by 
an inhomogeneous optical depth to the last scattering surface 
(Dvorkin & Smith 2009), 



^0 jv^fi ^^(2 )■ 



(14) 



Note that one might in principle be concerned about order 
corrections to the covariance matrix, particularly at high £ (see, 
e.g., Notari & Quartin 2012). However, these are small, provided 
that the spectra are relatively smooth. Although the order /3 (- 
10"^) deflections give large changes to the flf,„s for £ > 10^, the 
changes to the overall covariance are small (Chluba 201 1), since 
the deflection effect is coherent over very large scales. 

Note that these boosting-induced correlations couple £ with 
{ + I modes. They therefore share this property with pure dipo- 
lar amplitude modulations studied in the context of primordial 
statistical anisotropy (see, e.g, Prunet et al. 2005), as well as 
with dipolar modulations in more general physical parameters 
(Moss et al. 2011). However, these other cases do not share the 
frequency dependence of the boosting modulation effect, since 
they are not accompanied by a temperature dipole. 

We measure the peculiar velocity dipole using quadratic es- 
timators, essentially summing over the covariance matrix of the 
observed CMB fluctuations, with weights designed to optimally 
extract fi. A general quadratic estimator x for /3 is given by 
(Hanson & Lewis 2009) 



1 



{ 



xlm[T] = -Nf 2 Z Z 



M 



£1 £2 L 
nix ^2 ~M 



Ci=lmm h^lmm '«1.™2 

^{TtunTt2.n._-{f,^„T,^,„S), (15) 

where T(,„ are a set of inverse-variance filtered temperature mul- 
tipoles, W'^^g^^ is a weight function and N'^^' is a normalization. 
The ensemble average term () is taken over signal-i-noise real- 
izations of the CMB in the absence of velocity eff'ects. It cor- 
rects for the statistical anisotropy induced by effects like beam 
asymmetry, masking, and noise inhomogeneity. We evaluate this 
term using Monte Carlo simulations of the data, as discussed in 
Sect. 4. 

We use three different quadratic estimators to measure the ef- 
fect of boosting. The first, 0, simply adopts the weight function 
^titiL' ^'^'^ provides a minimum-variance estimator of the to- 
tal peculiar velocity effect. The two additional estimators, and 
f, isolate the aberration and modulation aspects of the peculiar 
velocity effect, respectively. This can be useful, as they are qual- 
itatively quite distinct effects, and suffer from different potential 
contaminants. The modulation effect, for example, is degenerate 



with a dipolar pattern of calibration errors on the sky, while the 
aberration effect is indistinguishable from a dipolar pattern of 
pointing errors. 

There is a subtlety in the construction of these estimators, 
due to the fact that the covariances described by and are 
not orthogonal. To truly isolate the aberration and modulation 
effects, we form orthogonalized weight matrices as 



W: 



Kt.L-Kc.LK'II^L and 



where the response function H is given by 



^ (2L+1) 



2 Z 



(16) 
(17) 



(18) 



with x,z - Py, (p, T. The construction of these estimators is analo- 
gous to the construction of "bias-hardened" estimators for CMB 
lensing (Namikawa et al. 2012). The spectra F( are diagonal ap- 
proximations to the inverse variance filter, which takes the sky 
map T T. We use the same inverse variance filter as that used 
for the baseline results in Planck Collaboration XVII (2013), and 
the approximate filter functions are also specified there. Note 
that our estimator is slightly different from that used in Planck 
Collaboration XVII (2013), due to the fact that we have orthog- 
onalized it with respect to r. 

The normalization N'^^'' can be approximated analytically as 



(19) 



This approximation does not account for masking. On a masked 
sky, with this normalization, we expect to find that (xlm) = 

fLM,skyPLM, whcrC 

/LM,.sky ^ jdn Yl^ih )M(n ■ n . (20) 

Here M(n ) is the sky mask used in our analysis. For the fidu- 
cial sky mask we use (plotted in Fig. 2, which leaves approx- 
imately 70% of the sky unmasked), rotating /z,M,sky onto our 
three basis vectors we find that /||. sky = 0.82, /l, sky - 0.17 
and /x, sky = -0.04. The large effective sky fraction for the jSy 
direction reflects the fact that the peaks of the expected velocity 
dipole are untouched by the mask, while the small values of /sky 
for the other components reflects that fact that the masking pro- 
cedure does not leak a large amount of the dipole signal in the 
)8|| direction into other modes. 



4. Data and simulations 

Given the frequency-dependent nature of the velocity effects we 
are searching for (at least for the t component), we will focus 
for the most part on estimates of yS obtained from individual fre- 
quency maps, although in Sect. 6 we will also discuss the analy- 
sis of component-separated maps obtained from combinations of 
the entire Planck frequency range. Our analysis procedure is es- 
sentially identical to that of Planck Collaboration XVII (2013), 
and so we only provide a brief review of it here. We use the 
143 and 217 GHz Planck maps, which contain the majority of 
the available CMB signal probed by Planck at the high mul- 
tipoles required to observe the velocity effects. The 143 GHz 
map has a noise level that is reasonably well approximated by 
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45/zKarcmin white noise, while the 217 GHz map has approx- 
imately twice as much power, with a level of 60;uKarcmin. 
The beam at 143 GHz is approximately 7' FWHM, while the 
217 GHz beam is 5' FWHM. This increased angular resolu- 
tion, as well as the larger size of the r-type velocity signal at 
higher frequency, makes 217 GHz slightly more powerful than 
143 GHz for detecting velocity effects (the HEX 100 GHz and LFI 
70 GHz channels would offer very little additional constraining 
power). At these noise levels, for 70% sky coverage we Fisher- 
forecast a 20% measurement of the component at 217 GHz 
(or, alternatively, a 5cr detection) or a 25% measurement of 
yS|l at 143 GHz, consistent with the estimates of Kosowsky & 
Kahniashvili (2011). As we will see, our actual statistical error 
bars determined from simulations agree well with these expec- 
tations. 

The Planck maps are generated at HEALPix (Gorski et al. 
2005)* A^side = 2048. In the process of mapmaking, time-domain 
observations are binned into pixels. This effectively generates a 
pointing error, given by the distance between the pixel centre to 
which each observation is assigned and the true pointing direc- 
tion at that time. The pixels at A^side - 2048 have a typical dimen- 
sion of 1.7'. As this is comparable to the size of the aberration 
effect we are looking for, this is a potential source of concern. 
However, as discussed in Planck Collaboration XVII (2013), the 
beam-convolved CMB is sufficiently smooth on these scales that 
it is well approximated as a gradient across each pixel, and the 
errors accordingly average down with the distribution of hits in 
each pixel. For the frequency maps that we use, the rms pix- 
elization error is on the order of 0.1', and not coherent over the 
large dipole scales which we are interested in, and so we neglect 
pixelization effects in our measurement. 

We will use several data combinations to measure p. The 
quadratic estimator of Eq. (15) has two input "legs," i.e., the {\mi 
and ^2'W2 terms. Starting from the 143 and 217 GHz maps, there 
are three distinct ways we may source these legs: (1) both legs 
use either the individually filtered 143 GHz or 217 GHz maps, 
which we refer to as 143 x 143 and 217 x 217, respectively; (2) 
we can use 143 GHz for one leg, and 217 GHz for the other, re- 
ferred to as 143 X 217; and (3) we can combine both 143 and 
217 GHz data in our inverse-variance filtering procedure into 
a single minimum-variance map, which is then fed into both 
legs of the quadratic estimator We refer to this final combina- 
tion schematically as "143-1-217." Combinations (2) and (3) mix 
143 and 217 GHz data. When constructing the weight function 
of Eq. (12) for these combinations we use an effective by - 2.5. 
To construct the f, which are the inputs for these quadratic esti- 
mators, we use the filtering described in Appendix A of Planck 
Collaboration XVII (2013), which optimally accounts for the 
Galactic and point source masking (although not for the inho- 
mogeneity of the instrumental noise level). This filter inverse- 
variance weights the CMB fluctuations, and also projects out the 
857 GHz Planck map as a dust template. 

To characterize our estimator and to compute the mean-field 
term of Eq. (15), we use a large set of Monte Carlo simula- 
tions. These are generated following the same procedure as those 
described in Planck Collaboration XVII (2013); they incorpo- 
rate the asymmetry of the instrumental beam, gravitational lens- 
ing effects, and realistic noise realizations from the FFP6 sim- 
ulation set described in Planck Collaboration I (2013). There 
is one missing aspect of these simulations which we discuss 
briefly here: due to an oversight in their preparation, the gravita- 
tional lensing component of our simulations only included lens- 
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ing power for lensing modes on scales L > 2, which leads to a 
slight underestimation of our simulation-based error bars for the 
(p component of the velocity estimator The lensing dipole power 
in the fiducial ACDM model is C^^ ^ 6 x 10"^, which repre- 
sents an additional source of noise for each mode of )8, given 
by cr^ iens = 1.2 X 10""^, or about one tenth the size of the ex- 
pected signal. The part of the estimator contributes approx- 
imately 46% of the total [i estimator weight at 143 GHz, and 
35% at 217 GHz. Our measurement errors without this lensing 
noise on an individual mode of p are erg ^ 2.5 x 10 while 
with lensing noise included we would expect this to increase to 

^(7-2 + (4/10)2o-2j^^^ = 2.54 X 10""^. This is small enough that 

we have neglected it for these results (rather than include it by 
hand). 

We generate simulations both with and without peculiar ve- 
locity effects, to determine the normalization of our estimator, 
which, as we will see, is reasonably consistent with the analyt- 
ical expectation discussed around Eq. (20). All of our main re- 
sults with frequency maps use 1000 simulations to determine the 
estimator mean field and variance, while the component separa- 
tion tests in Sect. 6 use 300 simulations. 

5. Results 

We introduce our results visually in Fig. 3, where we plot the 
total measured dipole direction )8 as a function of the maximum 
temperature multipole ^n,ax used as input to our quadratic estima- 
tors. We can see that all four of our 143/217 GHz based estima- 
tors converge toward the expected dipole direction at high ^^ax- 
At /"max < 100, we recover the significantly preferred direction 
of Hoftuft et al. (2009), which is identified when searching for a 
dipolar modulation of the CMB fluctuations. The t component 
of the velocity effect is degenerate with such a modulation (at 
least at fixed frequency), and (p gets little weight from i < 100, 
so this is an expected result. The significance of this preferred di- 
rection varies as a function of smoothing scale (Hanson & Lewis 
2009; Bennett et al. 2011; Planck Collaboration XXIII 2013). 
To minimize possible contamination on our results by this po- 
tential anomaly, from here onward we restrict the temperature 
multipoles used in our yS estimation to (^i^ - 500. This cut re- 
moves only about 10% of the total number of modes measured 
by Planck, and so does not significantly increase our error bars. 

There is a clear tendency in Fig. 3 for the measured velocity 
to point toward the expected direction jSy. At Planck noise levels, 
we expect a one standard-deviation uncertainty on each compo- 
nent of p of better than 25%. A 25% uncertainty corresponds to 
an arctan(l/4) = 14° constraint on the direction of p. We plot 
this contour, as well as the corresponding two standard devia- 
tion contour arctan(2/4) = 26°. It is apparent that the measured 
velocity directions are in reasonable agreement with the CMB 
dipole. 

We now proceed to break the measurement of Fig. 3 into 
its constituent parts for {j^va - 2000 (and truncating now at 
(mm - 500). In Fig. 4 we plot our quadratic estimates of the 
three components of p, as well as the decomposition into aberra- 
tion and modulation components for each of our four frequency 
combinations. The vertical lines in Fig. 4 give the amplitude es- 
timates for each component measured from the data, while the 
coloured and grey histograms give the distribution of these quan- 
tities for the 143x217 estimator for simulations with and without 
velocity effects, respectively (the other estimators are similar). 
As expected, the velocity effects show up primarily in ySy; there 
is little leakage into other components with our sky mask. For all 
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Fig. 3. Measured dipole direction p in Galactic coordinates as a function of the maximum temperature multipole used in the 
analysis, {max- We plot the results for the four data combinations discussed in Sect. 4: 143 x 143 (T symbol); 217 x 217 (a symbol); 
143 X 217 (x symbol); and 143 + 217 (+ symbol). The CMB dipole direction P\\ has been highUghted with 14° and 26° radius 
circles, which correspond roughly to our expected uncertainty on the dipole direction. The black cross in the lower hemisphere is 
the modulation dipole anomaly direction found for WMAP at f^ax = 64 in Hoftuft et al. (2009), and which is discussed further in 
Planck Collaboration XXIII (2013). Note that all four estimators are significantly correlated with one another, even the 143 x 143 
and 217 x 217 results, which are based on maps with independent noise realizations. This is because a significant portion of the 
dipole measurement uncertainty is from sample variance of the CMB fluctuations, which is common between channels. 



four estimators, we see that the presence of velocity along ySy is 
strongly preferred over the null hypothesis. At 143 GHz this sig- 
nal comes from both and fy. At 217 GHz it comes primarily 
from f||. Additionally, there is a somewhat unexpected signal at 
217 GHz in the j8x direction, again driven by the t component. 
Given the apparent frequency dependence, foreground contami- 
nation seems a possible candidate for this anomalous signal. We 
will discuss this possibility further in the next section. 

In Table 1 we present values for the j8 measurements of 
Fig. 4 under both the null hypothesis of no velocity effects, as 
well as assuming the expected velocity direction and amplitude. 
We can see that all of our measurements are in significant dis- 
agreement with the "no velocity" hypothesis. The probability- 
to-exceed (PTE) values for the "with velocity" case are much 
more reasonable. Under the velocity hypothesis, 217 x 217 has 
the lowest PTE of 1 1%, driven by the large 0x- 

In Table 2 we focus on our measurements of the velocity 
amplitude along the expected direction ySy, as well as perform- 
ing null tests among our collection of estimates. For this table, 
we have normalized the estimators, such that the average of ySy 
on boosted simulations is equal to the input value of 369 km s"' . 
For all four of our estimators, we find that this normalization 
factor is within 0.5% of that given by A^"'^*'/||,sky, as is already ap- 



parent from the triangles along the horizontal axis of Fig. 4. We 
can see here, as expected, that our estimators have a statistical 
uncertainty on ySy of between 20% and 25%. However, several 
of our null tests, obtained by taking the differences of pairs of /3\\ 
estimates, fail at the level of two or three standard deviations. We 
take the 143 x 217 GHz estimator as our fiducial measurement; 
as it involves the cross-correlation of two maps with indepen- 
dent noise realizations it should be robust to noise modelling. 
Null tests against the individual 143 and 217 GHz estimates are 
in tension at a level of two standard deviations for this estima- 
tor We take this tension as a measure of the systematic differ- 
ences between these two channels, and conservatively choose 
the largest discrepancy with the 143 x 217 GHz estimate, namely 
0.31, as our systematic error We therefore report a measurement 
of V|| = 384 kms"' + 78kms"' (stat.) + 115 kms"' (syst.), a sig- 
nificant confirmation of the expected velocity amplitude. 

6. Potential contaminants 

There are several potential sources of contamination for our es- 
timates above which we discuss briefly here, although we have 
not attempted an exhaustive study of potential contaminants for 
our estimator. 
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Fig. 4. Measurements of fi using combinations of the 143 and 217 GHz Planck maps, normalized using Eq. (19) and then divided 
by the fiducial amplitude of /3 - 1.23 x lO"-'. These estimates use {^[^ - 500 and /"max - 2000. In addition to the total minimum 
variance estimate 0, the measurement is also broken down into its aberration-type part, 0, and modulation-type part, f. Vertical lines 
give the Planck measurement for the four estimates described in the text. Grey histograms give the distribution of estimates for 
simulations of the 143 x 217 estimator, which do not contain peculiar velocity effects (the other estimators are very similar). The 
red histograms give the distribution for simulations which do contain peculiar velocity effects, simulated with the fiducial direction 
(along jS||) and amplitude. Black triangles on the x-axis indicate the relevant component of /sky given by Eq. (20), which agrees well 
with the peak of the velocity simulations. 
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Fig. 5. Plot of velocity amplitude estimates, similar to Fig. 4, but using an array of component-separated maps rather than specific 
combinations of frequency maps. The production and characterization of these component separated maps is presented in Planck 
Collaboration XII (2013). Histograms of simulation results without velocity effects are overplotted in grey for each method; they 
are all very similar. Vertical coloured bars correspond to the maps indicated in the legend, using the combination of our fiducial 
galaxy mask (which removes approximately 30% of the sky) as well as the specific mask produced for each component separation 
method. We see significant departures from the null-hypothesis simulations only in the jSy direction, as expected. Vertical black lines 
show the 143 x 217 measurement of Fig. 4. Please note the discussion about the subtleties in the normahzation of these estimates in 
Sect. 6. T 
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Table 1. Significance measures for the p estimates of Fig. 4. We 
form ?lx^ for the three measured modes of fi, using the mean and 
covariance matrix measured for simulations either with or with- 
out velocity effects. The covariance matrices are very similar in 
both cases, and so the difference between these two cases is only 
the mean-field subtraction. The "PTE" columns give the corre- 
sponding probability-to-exceed values for a distribution with 
three degrees of freedom. The measured values are extremely 
unlikely under the "no velocity" hypothesis, but very compatible 
for the "with velocity" case. 







P significance 




No velocity 


With velocity 






PTE [%] 


x' 


PTE [%] 


143 X 143 . . 


. 21. \ 


0.0005 


1.9 


58.87 


217x217 . . 


. 20.7 


0.0123 


6.0 


11.18 


143 + 217 . . 


. 25.1 


0.0015 


3.3 


35.44 


143x217 . . 


. 27.6 


0.0005 


1.8 


62.29 



Table 2. Measured velocity amplitudes along the jSy direction, in 
units of 369 km s \ using combinations of the 143 and 217 GHz 
data, as discussed in Sect. 4. The diagonal shows the results for 
the indicated reconstruction. Below the diagonal, the numbers 
given are for the difference of the two results, and the uncertainty 
accounts for the correlation between each pair of measurements. 
This lower triangle is a null test. Several of these null tests fail 
between 143 and 217 GHz, although both channels provide ev- 
idence for y6|| consistent with the expected amplitude and dis- 
crepant with zero at between four and five standard deviations. 







y8|| Amplitude 






143 X 143 


217x217 143 + 217 


143x217 


143 X 143 . . 
217x217 . . 
143 + 217 . . 
143x217 . . 


. 1.35 ±0.26 
. 0.60 ±0.21 
. 0.43 ±0.15 
. 0.31 ±0.13 


0.75 ±0.19 
-0.17 ±0.09 0.92 ±0.20 
-0.29 ±0.12 -0.12 ±0.07 


1.04 ±0.21 



Galactic Foregrounds: Given the simplicity of the fore- 
ground correction we have used (consisting only of masking the 
sky and projecting out the Planck 857 GHz map as a crude dust 
template), foreground contamination is a clear source of con- 
cern. The frequency dependence of the large fi^ signal seen at 
217 GHz, but not at 143 GHz, seems potentially indicative of 
foreground contamination, as the Galactic dust power is approx- 
imately 10 times larger at 217 than at 143 GHz. To test the pos- 
sible magnitude of residual foregrounds, we apply our veloc- 
ity estimators to the four component-separated CMB maps of 
Planck Collaboration XII (2013), i.e., NILC, SMICA, SEVEM, and 
RULER. Each of these methods combines the full set of nine 
Planck frequency maps from 30 to 857 GHz to obtain a best- 
estimate CMB map. To characterize the scatter and mean field of 
each method's map we use the set of common simulations which 
each method has been applied to. These simulations include the 
effect of the aberration part of the velocity dipole, although not 
the frequency-dependent modulation part. For this reason, it is 
difficult to accurately assess the normalization of our estimators 
when applied to these maps, particularly as they can mix 143 
and 217 GHz as a function of scale, and the modulation part is 
frequency dependent. We can, however, study them at a quali- 



tative level. The results of this analysis are shown in Fig. 5. To 
construct our p estimator for the component separated maps we 
have used by - 2.5, assuming that they contain roughly equal 
contributions from 143 and 217 GHz. Note that, because the sim- 
ulations used to determine the mean fields of the component sep- 
arated map included the aberration part of the velocity effect, 
it will be absorbed into the mean field if uncorrected. Because 
the aberration contribution is frequency independent (so there 
are no issues with how the different CMB channels are mixed), 
and given the good agreement between our analytical normaliza- 
tion and that measured using simulations for the frequency maps, 
when generating Fig. 5 we have subtracted the expected velocity 
contribution from the mean field analytically. We see generally 
good agreement with the 143 x 217 estimate on which we have 
based our measurement of the previous section; there are no ob- 
vious discrepancies with our measurements in the ySy direction, 
although there is a somewhat large scatter between methods for 
011 . In the )8x direction the component-separated map estimates 
agree well with the 143 x 217 estimator, and do not show the 
significant power seen for 217 x 217, suggesting that the large 
power which we see there may indeed be foreground in origin. 

Calibration errors: Position-dependent calibration errors in 
our sky maps are completely degenerate with modulation-type 
effects (at fixed frequency), and so are very worrisome as a 
potential systematic effect. We note that the Planck scan strat- 
egy strongly suppresses map calibration errors with large-scale 
structure (such as a dipole). As the satellite spins, the detectors 
mounted in the focal plane inscribe circles on the sky with open- 
ing angles of between 83° and 85° (for the 143 and 217 GHz de- 
tectors we use). For a time-dependent calibration error to project 
cleanly into a dipolar structure on the sky, it would need to 
have a periodicity comparable to the spin frequency of the satel- 
lite (Imin"'). Slower fluctuations in the calibration should be 
strongly suppressed in the maps. There are calibration errors 
which are not suppressed by the scan strategy, however. For ex- 
ample, a nonlinear detector response could couple directly to 
the large CMB dipole temperature. Ultimately, because position- 
dependent calibration errors are completely degenerate with the 
T component of the velocity effect, the only handle which we 
have on them for this study is the consistency between and 
f . From another viewpoint, the consistency of our measurement 
with the expected velocity modulation provides an upper bound 
on dipolar calibration errors. 

Pointing errors: In principle, errors in the pointing are per- 
fectly degenerate with the aberration-type velocity effect in the 
observed CMB. However the Planck pointing solution has an un- 
certainty of a few arcseconds rms in both the co- and cross-scan 
directions (Planck Collaboration VI 2013). The 3' rms aberra- 
tions induced by velocity effects are simply too large to be con- 
taminated by any reasonable pointing errors. 

7. Conclusions 

From Fig. 3 it is clear that small-scale CMB fluctuations ob- 
served by the Planck data show evidence for velocity effects in 
the expected direction. This is put on more quantitative footing 
in Fig. 4 and Table 2, where we see that all four of the 143 and 
217 GHz velocity estimators which we have considered show ev- 
idence for velocity effects along ySy at the level of more than four 
standard deviations. Detailed comparison of 143 and 217 GHz 
data shows some discrepancies, which we have taken as part of 
a systematic error budget; however, tests with component sepa- 
rated maps shown in Fig. 5 provide a strong indication that our 
217 GHz map has slight residual foreground contamination. The 



8 



Planck Collaboration: Doppler boosting of the CMB: Eppur si muove 



component separated results are completely consistent with the 
143 X 217 estimator which we quote for our fiducial result. 

Beyond our peculiar velocity's effect on the CMB, there 
have been many studies of related effects at other wavelengths 
(e.g., Blake & Wall 2002; Titov et al. 2011; Gibelyou & Huterer 
2012). Closely connected are observational studies which ex- 
amine the convergence of the clustering dipole (e.g., Itoh et al. 
2010). Indications of non-convergence might be evidence for 
a super-Hubble isocurvature mode, which can generate a "tilt" 
between matter and radiation (Turner 1991), leading to an ex- 
tremely large-scale bulk flow. Such a long-wavelength isocur- 
vature mode could also contribute a significant "intrinsic" com- 
ponent to our observed temperature dipole (Langlois & Piran 
1996). However, a peculiar velocity dipole is expected at the 
level of yfC\ ~ 10"^ due to structure in standard ACDM (see, 
e.g., Zibin & Scott 2008), which suggests a subdominant intrin- 
sic component, if any. In addition, such a bulk flow has been sig- 
nificantly constrained by Planck studies of the kinetic Sunyaev- 
Zeldovich effect (Planck Collaboration Int. XIII 2013). In this 
light, the observation of aberration at the expected level reported 
in this paper is fully consistent with the standard, adiabatic, pic- 
ture of the Universe. 

The Copernican revolution taught us to see the Earth as or- 
biting a stationary Sun. That picture was eventually refined to 
include Galactic and cosmological motions of the Solar System. 
Because of the technical challenges, one may have thought it 
very unlikely to be able to measure (or perhaps even to define) 
the cosmological motion of the Solar System. . . and yet it moves. 
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